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ABSTRACT 


Th« problem of roll lock-in has gained importance with the 
use of finned projectiles as anti-tank ammunition. We 
consider a slightly asymmetric projectile with a trim angle 
that rotates with the projectile and shows large 
amplification when the projectile roll rate resonates with 
its nutation frequency. For a projectile with an offset 
cent er-of-mass , the trim angle induces a roll moment which 
can cause the roll rate to lock in at resonance. 

The present study considers lock-in solutions obtained for 
different combinations of parameters of the nonlinear system 
for the projectile angular motion. Both normal lock-in, 
where the spin is in the same sense as the expected steady- 
state spin, and reverse lock-in, where the spin is in the 
opposite sense, are shown. Stability of these lock-in 
solutions and the influence of initial conditions on the 
steady-state response are studied. Stable lock-in is seen to 
occur only for a small fraction of the initial conditions 
tested. This is shown to be due to the presence of an 
unstable equilibrium near a stable lock-in solution. 

Other anomalous roll behavior of rolling finned projectiles 
like transient resonance, catastrophic yaw, and roll break- 
out are also studied. The present analysis reveals that a 
projectile at transient resonance shows increased damping in 
yaw, contrary to the expected yaw build-up due to resonance. 
It is shown that catastrophic yaw can be explained without 
introducing nonlinear side moments in the analysis. 
Furthermore, it is shown that lock-in can occur even when 
the equilibrium at resonance is known to be unstable. 
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Chapter 1 


INTRODUCTION 

Finnad projectiles have today gained common use as military 
ordnance, especially as anti-tank ammunition. Such 
projectiles (high Kinetic Energy weapons) are fired at very 
high velocities (over 1500 m/s) and are intended to pierce 
tank armour by their sheer speed, and consequent energy. 
Given the need for high muzzle velocities, these projectiles 
are preferred to be fin-stabilized in contrast to spin- 
stabilized artillery shells and projectiles. The use of 
fin-stabilization has, however, brought into focus a number 
of problems which, although identified long ago, were not 
considered of sufficient practical significance. Among these 
are roll resonance and the associated problem of roll locfc- 
i n . 

Finned projectiles are. designed to be inherently statically 
stable and, hence, do not need to be spun. Such an unspun 
projectile which is symmetric exhibits precessional and 
nutational motion. Host projectiles of this type are 
dynamically stable, that is, the precessional and nutational 
modes damp out with time. However, slight configurational 
asymmetries always exist . "Slight configurational asymaetry" 
refers to missiles where the assumption of rotational 
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symmetry is valid for JLlC-jd e r-i.v.a t i v e s a n d moments of 
inertia; however, significant asymmetric aerodynamic moments 
are generated so as to trim the unspun projectile at an 
angle of attack (sometimes called the "asymmetry angle"). 

Projectiles are usually designed with marginal static 
stability that results in small natural frequencies in pitch 
and yaw. Moreover, small asymmetries or disturbances can 
easily set an initially unspun projectile into a roll. Such 
an accidentally produced roll rate may readily match with 
the low natural pitch-yaw frequency. This results In 
magnification of the trim angle, and the projectile 
undergoes large-amplitude yawing motion. This is called roll 
resonance and may lead to instability. 

A common solution to avoid roll resonance is to cant the 
fins and deliberately roll the projectile through the 
resonance region. Surprisingly, many flight failures were 
observed when the roll rate failed to build up to the 
equilibrium (design) value corresponding to the fin cant. 
Instead, the roll rate was seen to get locked in at 
resonance; this phenomenon is termed roll lock-in. The 
observed lock-in tendencies could not be explained based on 
the existing linear aerobailistic theory. The theory 
predicted that the roll rate would vary through the 
resonance region and the projectile would attain only a 
fraction of the maximum (resonant) amplitude. The actual 
disturbances were much larger. To satisfactorily explain the 
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occurrence of roll locK-in and to account for the 
discrepancies between the linear theory and actual flight 
experiences, various nonlinear formulations have been 
suggested in the literature. The present study is also a 
step in this direction. 

Before proceeding with a summary of the scope and extent of 
the present study, a brief survey of existing literature is 
presented. A fairly complete linear theory of the pitching 
and yawing motion of spinning shells was first published by 
Fowler et alCtl in 1920. This was later refined by McShane 
et al C£1 in 1953. The linear theory has been used in this 
form to date and a good description is available in 
textbooks [13, 41 . A recent survey paper by Murphy ESI serves as 
an excellent guide to the field of “flight dynamics of 
missiles." It concentrates on various sources of dynamic 
instability and is of direct relevance to the present study. 
A survey paper by PlatusC61 on ballistic re-entry vehicle 
dynamics also covers many related topics. 

The concept of a "slightly asymmetric missile” was 
introduced by NicolaidesCTl in 1953. He modified the linear 
theory and applied it to analyze angular motion of such 
missiles and showed the possibility of occurrence of roll 
resonance. Murphy C81 predicted the influence of nonlinear 
forces and moments on a symmetric missile. The phenomena of 
"roll lock-in” and "catastrophic yaw" were put forward by 
Nicolaides t93 soon after. The mechanism of generation of 
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nonlinear noments causing these phenomena was explored by 
PriceClOD, ChadwickElIl, and BarberaCIEl. However, most of 
the effort was directed towards re-entry vehicles. 

A more general analysis of the motion of asymmetric rolling 
bodies in the presence of nonlinear aerodynamics was 
presented by Murphytl33. Using the method of slowly-varying 
parameters, ClareC143 was able to extend his solution to 
include resonance. However, both t1urphyCt33 and ClareC143 
had assumed a constant spin rate. This constraint was later 
relaxed by Murphy[153, but the rolling moment was assumed to 
be linear. Thus, none of these efforts could explain roll 
lock-in. 

Nayfeh and SaricC1£3 employed the method of multiple scales 
and solved the coupled nonlinear roll, pitch, and yaw 
equations to analyze roll resonance of a re-entry vehicle 
and to obtain necessary conditions for lock-in. Modeling and 
analysis of roll resonance and lock-in of re-entry vehicles 
was also taken up by KevorkianEI 73 , and Kevorkian and 
LewinE183. However, a recent work by Murphy Cl 93 has gone a 
long way in presenting a lucid exposition of the theory of 
roll lock-in of finned projectiles. 

The present study seeks to expand on the work of Murphy Cl 93 . 
The question of stability of the equilibrium solution at 
lock-in obtained from the nonlinear system equations is 
considered in greater detail . Lock-in solutions for the case 
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of an offset cent re-of-mass are investigated for varying 
sets of parameters. The importance of both stable and 
unstable lock-in solutions is brought out. Other anomalous 
roll behavior of rolling finned projectiles like transient 
resonance is also considered. Finally, we take a look at 
the phenomenon of catastrophic yaw and at the possibility of 
break-out from lock-in of finned projectiles. 

The remainder of this thesis is arranged as follows. Chapter 
£ se^s down the equations of motion of a rolling projectile 
in an aeroba 1 1 i st i c coordinate frame. This is followed by a 
brief description of the main results of the linear 
aerobal 1 ist ic theory. Chapter 3 discusses the lock-in theory 
and investigates the stability of lock-in solutions. The 
results of the present study covering the phenomena of roll 
resonance and lock-in, transient resonance, catastrophic 
yaw, and roll break-out have been put together in Chapter 4. 
The thesis ends with a summary of its chief contributions, 
and with suggestions for future work based on these 
f oundat ions . 



Chapter II 


LINEAR THEORY 


£.1 Aerobal 1 ist ic Axes 

Aerobal 1 i St ic axes (Fig.S.t) provide a convenient non- 
spinning coordinate systen for the study of rolling 
missiles. The X-axis faces forward along the projectile 
principal axis while the Y-axis lies perpendicular to the 
X-axis and is initially in the horizontal plane. The'z'-axis 
is at right angles to both the X and Y axes, and points 
initially in the downward direction so as to for* a right- 
handed coordinate system. The axis system pitches and yaws 
with the projectile, but has a zero roll rate. Then, after 
the initial instant, the 'v-axis may be displaced from the 
horizontal plane by the pitching and yawing motion of the 
pro j ect i le C51 . The missile-fixed axis system (X,Y,Z) is also 
shown inFig.2.1. 

2.2 Equations of Motion 

The equations of motion for a basically symmetric 
missileC191 in the aerobal listic system are as follows: 

F^ = m ( ii + wq - rv ) 

Fy * m ( V + ru ) 


■ m ( w - qu ) 


. . . ( 2 . 1 ) 



xz 

M 1 SSILE -F IXED 

AXIS 

YZ 

AERoBALL 1 ST 1 C 

A KiS 



AFT, looking forward 


FIG. 2-1: AEROBALLISTIC AXIS SYSTEM 




B 


L=Ij^p ...<£.£) 

M »= I q + Ixpr 

N * I f - Ijcqp ...{£. 3) 

Linear aerobal 1 ist ic theory seeks to solve the equations of 
angular motion of a projectile under the assumptions of 
£i) constant roll rate, 

(ii) small angle of attack and linear aerodynamics, 

( i i i ) const ant total projectile velocity. 

In addition, the force equations (£.1) are taken to be 
identically satisfied. 

For a symmetric missile, 

• • 

q » r » * ...icL.4) 

where o( » w/u and ^ = v/u. 

The following expansions for the pitching and yawing moments 
are assumed C3,43 : 


"i"- V * Vi* 


In addition, a body-fixed asymmetry will cause a trim moment 
Hy which will rotate with the projectile roll rate. <£ 3) 
will now appear as under: 

Iq + rplx - - M^q - pj? = M.j.cos(pt ) 

if - qplx - - - N^r - N^e<pS * llj.sint pt )...(£. 6) 
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A complex angle of attack is defined as ^ *(2>+ i wi . 


Also, 

t he 

f ol lowing 

equal it ies 

hold 

due 

to rotational 

symmetry : 









«o( 



Nr 







- -N * 






(£.71 

<£. 6) 

can now be combined by 

multiplying 

the 

second 

equat ion 

by ' i 

* and 

adding it 

to the 

first 

to give a 

single 

equat ion 

in the 

complex variable^. 







• » 


-+ 

H^c)f 

- (- 


_ 



s * 

V r j 

X 


\ 

X 

xy 


= C JjL -cJcp ( L pi) 

^ ...( 2 . 8 ) 


Since p has been assumed constant, (£.6) is a linear 
differential equation with time-invariant coefficients and 
can be solved for the complex angle of attack,^. The trim 
moment appears on the right hand side of (£.81 and acts as a 
source of excitation. 

£.3 Tricyclic Motion 
Rewriting (£.8) as 

5 + N,| + N_2'i» N3exp(ipt} ...(£.9) 

a solution can be found as t33 

g » K.,exp(a^ t) + K^exp(m 2 t) + K^exp (ipt )..{£. 10 1 



where 


1 0 


m 


1.2 



+ 



-2. 


N. 


« >v 


1.2 




. . 1 1 J 


Kj and are obtained from the initial conditions while 
is given byC33 

^ OP~‘^OC 




Ns 

# O f o \ 


[_l( 


. . 1 c . 1 c 1 

(£.11) on 

subst itut ing 

for and fron (£.8) 

gives 

'^1,2 

/ lM<v + 

) ( 1 + a) -+• ^ 

. . (£. 13) 

” V 2 X ' 

Xa 


- 

* 

. . (£.14) 


2 X ^ 



where A 

J'-i. 

and s.* 

^ ( ^o</x) 

. . (E. 15) 


For a statically stable (finned) project ile ; ^ 

s^< 0, 0 < A < 1 - 

Therefore/ uj, is in the same sense as p and is larger in 
magnitude than uj^ which is in a sense opposite to p. uj, is 
called the nutation frequency and UJ 2 _ the precession 
frequency. The angular notion described by (£.10) consists 
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of three modeE and is represented by three rotating vectors 

(Fig.e.e) : 

NUTATION vector : ■ K,q expC(A, + iuj,)tJ 

"= to,t + arg ) 

PRECESSION vector: » K 2 Q*xp[(A_^+ ico^)t3 

0 ■ uj„ t + arg ((T ) 

TRIM vector : » K^^explipt) 

9^3 “ 

. . Ce. 16 ) 

The locus of the projectile tip in Fig.E.S is sometines 
called a yaw rosette. Examples of yaw rosettes are available 
in Ref C43 . 


£.4 Resonance 

Since rotates in a direction counter to the projectile 
roll motion, resonance is possible only for pe , as seen 
from (E.IE). For p “ uj, ,(E.14) yields: 



which gives the resonance roll rate as 

p « — ^ . . (E. 17) 






t3 


Since l^/l << 1 usually, the resonance roll rate can be seen 
to be very nearly equal to the projectile natural pitching 


frequency, - 

M^/I . When the 

project i le 

roll rate 

is 

near 

the resonance 

frequency , 

the yaw 

rosette 

has 

an 

approx imately 

circular pattern. Angular 

mot 1 on 

of 

this 


nature has been called lunar yawing notionC43. 

Projectiles at resonance show bounded oscillations of large 
amplitude due to an Increase in the magnitude of the trim 
angle of attack Using (2.12), the amplification of 

the non-rolling trim can be expressed asCl^l 

..(2.18) 


. . ( 2 . 19 ) 


- Mc< / J 

SE — 

*<3 I f t 

‘ (^ *■0 

where hj * ^2^ P - (jo^ ) + p ( ujj - pj 
h^, * \( p ) + A, { - p ) 


The magnitude of the trim amplification is obtained from 
(2.18). 




. . ( 2 . 20 ) 


This can be seen from (2.19) and (2.20) to grow to a maximum 
at resonance (p = ). Fig. 2. 3 taken from RefCSl shows the 
trim amplification at resonance. Thus, the linear theory is 
adequate to explain the occurrence of roll resonance. 




ol 0-5 1-0 VS 


Pig 2.3= TRIM AMPLIFICATION 
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In predicting the behavior of a projectile rolling through 
the resonance region, the variation of the roll rate cannot 
be overlooked. The linear theory, on the other hand, assuineE 
a constant roll rate without which (£.8) will not be 
amenable to an analytic solution in a closed form. However, 
it seems plausible to apply the linear theory to the problem 
of passage through resonance of a rolling projectile by 
assuming the roll rate to be varying slowly. Then, one can 
use the roll rate at every instant as a constant in the 
linear pitch-yaw equat i on (£ . 6 ) , and thereby obtain the 
resulting angle of attack at each instant given by <£.10). 
This approach, however, neglects any lag in the response due 
to the presence of damping forces. Notwithstanding this 
drawback, the above approach appears adequate to predict the 
behavior of a projectile in flight as it passes through 
resonance. 

But as flight experiences showed £1 0 , 1 1 3 , predictions based 
on the linear theory with the above modification were proven 
to be unsatisfactory. The theory was, in many cases, unable 
to explain observed flight behavior like roll lock-in and 
catastrophic yaw. This gives rise to the need to develop a 
more comprehensive theory. To that purpose, the assumption 
of constant roll rate is relaxed and some nonlinear terms 
are introduced in the linear equations of motion. This 
proves effective to explain the occurrence of roll lock-in 
as discussed in the following chapter. 



Chapter III 


LOCK-IN THEORY 

3.1 Roll Equat ion 

The main shortcoming of the linear theory has been seen to 

be its assumption of a constant roll rate. This must, 

therefore, be relaxed as the first step in the development 
of a lock-in theory. That requires the rolling moment in 
(£.E) to be non-zero. The rolling moment in its simplest 
form can be put down as comprising of two components: a 
roll-damping term, and a term due to fin cant. 

L « L^^p + 6^ .. .<3.t> 

The pitch-yaw equat i on ( £ . 8 ) retains its original form while 
the roll equation (£.£) with the rolling moment given by 

(3.1) is independent of g. Thus, it can be solved to give a 

steady-state roll rate 

Pss- “ tLs /L^) ...t3.2> 

Now, (£.8) can be solved as before with p»Pss Thus, only 
one equilibrium solution is possible. This is called the 
design solution. The value of Pg^iand the resultant ^1 can 
be selected by the designer by using an appropriate degree 
of fin cant, . 
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It may be observed that the expansion assumed in (3.1) for 
The rolling moment is insufficient to predict the occurrence 
of lock-in. An induced roll moment, which is a function of 
the total angle of attack and the roll angle between a 
particular fin and the plane of the total angle of attack, 
has been suggested [93 to explain the phenomenon of roll 
lock-in. With this added term, (3.1) modifies to 

L = Lj^p + Lg(fl,J) . . .(3.3) 

where Lg is the induced roll moment term. Now, (£.£) and 
{£.3) constitute a set of nonlinear equations and, in 
general, can have more than one steady-state equilibrium 
solution. 

3.£ Center-of-Mass Offset 

One source of an induced roll moment is a radially offset 
center-of-mass . Fig. 3.1 shows the center-of-mass offset by a 
quantity (lij^) in the <|>“0 plane (along the body-fixed Y- 
axis). The resultant rolling moment is given by 

Lg = Fj^. (if^^ ).sin • . . . C3.4> 

The differential equations for the angular motion of a 
projectile can now be derived. A brief outline of their 
derivation is provided in Appendix A. Ue follow the 
notation of Ref C193. 

The roll equation obtained from (£.£) and (3.3) is 




FlG-3.1 : PROJECTILE WITH AN OFFSET 


CENTER -OF- MASS 




. . . (3.5) 


1 9 


0 - (jg- iG - K «» 0 


The complete equation for the pitching and yawing motion 
derived from (£.1) and (£.3) is 


CH + i (£-o-)03 ^ + ( 1-0-) Cl-f>^+i0h+ii>3 
« ( 1-ff1 jhfexpf i^^) 


. (3.6) 


The independent variable in equations (3.5) and (3.5) is a 
dimensionless time, T , defined in equation (A. IE) of 
Appendix A. By virtue of this definition, the roll rate in 
the above equations is scaled by the value of the resonance 
frequency. Thus, resonance is signified by a value of 
101^1. Similarly, the complex angle of attack, '5 , is 
obtained by scaling ^ by the value of its magnitude at 
resonance as shown in equation (A. 17) of Appendix A. The 
final term within the brackets in (3.5) stands for the 
induced roll moment. The excitation due to the trim in the 
pitch-yaw equation appears on the right hand aide of (3.6). 
The dot (.) is used here, and henceforward, to denote 
derivative with respect to T . 

3.3 Equilibrium Solutions 

The equilibrium solutions of (3.5) and (3.6) can be obtained 
by solving the following equations: 

(ig - (ii * 16 < 


. . . (3.7) 



eo 



lh| (C 

• • 

I — ($£ -t- c h 


. . . (3.8) 


« • % ' # * 

(3.7) and (3.8) are obtained by setting and ^ in (3.5) 
and (3.6) to zero. In the above equations, ‘h' represents 
the aerodynamic asymmetry with orientation and '6' 
stands for the magnitude of the center-of-mass offset. 


(3.7) and (3.8) draw attention to some very significant 
points . 


(i).If no configurational asymmetries exist, then the 
projectile will have a zero trim angle. This implies 
that h is zero. Therefore, there will be no 
amplification of trim at resonance and the complex angle 
of attack will have zero magnitude at equilibrium as 
seen from (3.8). It follows from (3.7) that the 
equilibrium roll rate is equal to the design value which 
implies that the problem of roll lock-in is ruled out. 


(ii).No center-of-mass offset means 6 is zero, which 

• • 

from (3.7) implies that 9^ is always equal to <8^ (the 
design roll rate). For a non-zero trim th+O), reBonance 
({4)1 = 1) may occur with resultant trim amplification 
given by (S.SO). But, the rolling moment in this case 
takes the form (3.1) and a lock-in solution requiring 
({>^.^.+1 is impossible. Physically, this implies that no 
s t eady-st at e equi librium solut i on e x ist s in the vi ci nit y 


of resonance . 
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<iii).Even when h end G are both non-zero, (3.7) and 
(3.8) need not have a lock-in solution for all 
combinations of values of the parameters h, 0^^ t «nd G. 

The above statements are illustrated by an example. We 
choose 

fs « 3.0, h « 0.1. 

Four values of 0^,) are considered. 

O", 90‘ , leo' , 270" . 

G is varied from 0 to 5.0 in steps of 0.1. To solve (3.7) 
and (3.8) for 0^, the substitution of (3.8) in C3.7) leads 
to a fifth-order algebraic equation in The real roots 

representing the equilibrium roll rates (^^l obtained by 
solving this equation for the range of values of G are shown 
in Fig. 3. 2. 

The following inferences can be drawn from the figure: 

A design solution, S, exists near the design value 
(^g=3.0) in each of the four cases for all values of G. 

s' 

Lock-in solutions (^^"=+1 ) exist only for values of C 
greater than a particular minimum value. This ainiaua 
value is different in each case. 

Lock-in solutions are of two types: normal resonance, N, 

» m 

(^g®. + 1) and reverse resonance, R, (♦^■’.-1). 

Both the N and R types of lock-in solutions always 


appear in pairs. 
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If an equilibrium solution exists at resonance# it may or 
may not be stable. The question of stability of equilibrium 
solutions is addressed next. 

3.4 Lock-in Stability 

We make use of Lyapunov's first method to investigate the 
stability of the equilibrium solutions of (3.5) and (3.6). 
Following Ref 1193 # we assume a small perturbation of an 
equilibrium state. 

$ * ■'^ 2 * ^^3 

■n4 

TJg. ...(3.9) 

(3.5) and (3.6) can be cast as a system of five first-order 
differential equations. In vector notation# 

i. » X(x) . . (3. t 0) 

Using (3.9)# (3.10) can be written as 

1^ » A ‘Tj + g(-q) ..(3.11)^^^^^^^^^^^^^^^^ 

where g_(5^) represents higher-order terms. The elements of 
the 5x5 matrix A# called the stability matrix# are listed in 
Table 3 . 1 Cl 93 . Incidental ly, some of the terms appearing in 
the stability matrix of RefC193 are erroneous and they have 
been corrected in Table 3.1. 



TABLE 3.1 STABILITY MATRIX 
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(3.t1) can ba approximated by its linear part to study its 
stability characteristics provided g(7)) is negligible, which 
is usually the case. In the linear approximation, 

m 

. . c3. ts) 

The condition for stability requires that all the 

eigenvalues of matrix A have negative real parts. 

For an investigation of equilibrium solutions and their 

stability, one value of G is selected for each of the four 
cases depicted in Fig.3.&. In addition to the earlier 
parameters, the following parameters are selected to have 
the values given below: 

Equilibrium solutions are obtained by solving (3.7) and 

(3.8) for the roll rate and complex angle of attack 
respectively. These are listed in Table 3.£ for all the four 
cases considered. The stability matrix, A, for an 

equilibrium state is obtained by evaluating the elements in 
Table 3.1. The stability of the equilibrium solution is, 
then, determined by the eigenvalues of the resulting matrix. 
The equilibrium solutions tabulated in Table 3.S are 
classified as under: 

SD Stable - Design solution 

SN : Stable - Normal Resonance solution 

SR : Stable - Reverse Resonance solution 


U 


Un s t a b 1 e s o 1 ut i o n 




G 

SOLUTIONS 

1 

CLASSIFICATION 


Ret rmfSe) 



0. 1, H « 0.1, Kj. 


PARAMETERS: » 3.0, h 


01 , 0 - - 0.1 . 
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The following observations can be made from Table 3.E: 

For every stable lock-in solutioni there exists an 
unstable solution whose equilibrium roll rate is very 
near that of the stable lock-in solution. 

The design solution in each case is stable. Thus, all 
the cases under consideration show two stable solutions. 
For a different set of parameters, the number of stable 
solutions obtained may be different from two. 

A nonlinear system may approach any one of its many possible 
stable states or diverge to infinity. Thus, for a nonlinear 
system with two stable solutions, it is impossible to 
predict a priori which of the two solutions will occur once 
the system is set in operation. In the context of the 
examples in Table 3.£, the projectile may attain the design 
solution or it may find itself at a stable resonance 
solution. This distinction is Influenced solely by one 
factor, viz. the initial conditions at launch. This is taken 
up in the next chapter where the conditions under which a 
projectile does get locked in are considered. 



Chapter IV 


RESULTS AND DISCUSSION 

In nonlinear eyatems, which one of the Many possible 
equilibrium solutions actually occurs depends on the initial 
conditions. The final equilibrium state attained for a 
certain set of initial conditions can be deternined only by 
studying the response of the system obtained by solving its 
equations of motion. For a given a set of initial 
conditions^ the response of a projectile in flight can be 
obtained from a numerical integration of the system of 
equations (3.10). The numerical scheme used in the present 
study is a fourth-order Runge-Kutta-Ci 1 1 method. Using the 
above procedure, we first investigate the occurrence of roll 
lock-in. Following this, the phenomena of transient 
resonance, catastrophic yaw, and roll break-out are studied 
for finned projectiles. 

4.1 Effect of Initial Conditions 

The initial conditions for a projectile are the roll rate, 
the complex angle of attack, and the complex angular 
velocity at launch. To study responses corresponding to all 
possible initial conditions is clearly impracticable. 
Instead, a realistic set of values was considered. The 



initial roll rate and the initial complex angle (5c>) 

were taken to be zeroCi93, while the Initial complex angular 
velocity was varied. 

4 • - 0 

\ •= iSpI exp(i0*) ...{4.1) 

where ©* is the orientation of the initial total angular 
velocity ( \ ) . 

The following values were chosen for | | : 

I « 0.1 , O.e, 0.5, 1 . 0 . 

For each of the above values, O'^was set at 
©* - o', 90*, 180 *, E70*. 

Thus, a total of sixteen different Initial conditions were 
tested for each case listed in Table 3.£. 

The results of the study are as given below. Cases (a), (b), 
(c), and (d) refer to the four cases in Table 3.2 
corresponding to Fig. 3. 2. 

In cases (a) and (c)| the design solution occurs for 
each of the sixteen assumed set of initial conditions. 

In case (b), the normal resonance solution was seen only 
for {111 other initial conditions led to 

the design solution. For coaparislon, two sets of 
Initial conditions were chosen, one leading to resonance 
.nd th. oth.r yielding the de.ign .elution. The value of 
|4„| ua. kept fi.ed at l.«. For O-O", »e get th. nor.al 



lock“ir» solution snd 8**1B0 givos the design solution. 
The roll (spin) response end the yaw response 
corresponding to tbwse initial conditions are shown in 
Fig. 4.1. The figure highlights two significant results. 
First, the totally different nature of the response for 
two different sets of initial conditions Cparaseters 
held constant) is brought out very clearly. Second, the 
yaw response shows the large Magnitude of the total 
angle of attack at equilibrium for the resonance 
solution. This is in contrast to the yaw response for 
the design solution which can be seen to damp out with 
time. 

In case (d), |5el-t.O. ©-E7o", leads to the reverse 

resonance solution while the design equilibrium state is 
reached for the rest of the initial conditions. The roll 
and yaw responses in this case for ® ***• shown in 

Fig. 4. 8. Different initial conditions for reverse lock- 
in ( 0 - 270 ) and for the design solution (0-90) in Fig. 4.2 
again indicate the vastly different nature of the 
response- As in case tb), similar inferences about 
nature of the response can be drawn from Fig. 4.2. 

Based on observations such as mentioned for cases 

(d) above, designers have been cautionedCSl to avoid 

resonance. Catastrophic effects of lot'ge angle of 
amplification, range shortening fro. incraased drag, and 
structural failure have been mentionedtei as the disastrous 
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eff«ct* of roll lock-in. Howevor, our study suggests that 
the occurrence of stable resonance solutions is not very 
common. For example, of the sixty-four sets of Initial 
conditions considered in the present study, only two lead to 
a lock-in solution. 

In order to explain why lock-in solutions occur so rarely, 
we recall that resonance solutions never occur in isolation 
and that an unstable solution always exists near every 
stable lock-in solution. While it is not possible to provide 
a topological portrait of the trajectories for a fifth-order 
system, many of the points can be made by using a two- 
dimensional phase-plane picture. As an example, consider 
Fig. 4. 3 which shows an unstable saddle point, P, near a 
stable focus, Q. The point Q can be taken to represent a 
stable resonance solution with an adjacent unstable solution 
shown by P. Along which trajectory a system will find Itself 
depends on the initial conditions. Clearly, only those 
trajectories that lie between A and A' can reach the stable 
focus. The closer P and 0 are to one another, the smaller 
will be the gap between A and A', and fewer trajectorieB 
will find their way to the stable focus, Q, This underlines 

the Importance of the nature and location of all solutions, 

stable and unstable, in a nonlinear analysis. 




FIG.4-3 : PHASE -PLANE REPRESENTATION 
OF LOCK-IN SOLUTION 
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4.2 Tr«n*i»nt R*»or>«nc* 

As concludtd from th* previous discussion, a projsctils nay 
lock in at resonance for only a few initial conditions. 
However, its notion may still be significantly affected as 
it passes through the resonance region and the effect 
depends on the time spent at resonance condition. Depending 
on the orientation of the asymmetries, the nonlinear roll 
moment in (3.4) may either speed up or slow down the rolling 
motion of a projectile passing through resonance. As a 
result, the roll rate could linger at the resonance value 
for some tine before building up to the design value. This 
phenomenon is termed transient resonance CtSl . In terms of 
Fig. 4. 3, this is represented by trajectories like B which 
pass near the resonance solution, Q, but are not attracted 
to it. It is conjectured that the projectile would 
experience yawing motion of increased amplitude during the 
period it is at resonance. 

To investigate transient resonance, the following values of 
the parameters are used: 

4 - 3 0* h - 0.1, 

Case A: C * 6 Case B* 6 « 3.0 
For case A, Q is zero and only one solution is possible (the 
cJesign solut ion) . For case B, B is 3.0. The stable solutions 
for this case are shown in Table 4.1. The eigenvalues 
obtained from the stability matrix of Table 3.1 are listed 



lAikl di 


HAfiU Lfl£K-ltl mUTlONS 


«>K 

5 

SOLUTIONS 

EIGENVALUES 



4 

R*( Sg) 

I m ( Se J 


90 

0 

3.000 

0.0000 

0.0000 

-0 . 042 ± i 1.69 

-0.057 + 1 3.80 

-0.100 


3.0 

I . 020 

-0.8474 

0 . 3357 

-0.085 ± 1 0.69 

-0.045 + i 1.91 

-0.038 



2.920 

-0.0005 

- 

0.0133 

-0.046 ii i 1 .81 

-0.056 + i 3.73 

-0.093 

870 I 

0 

3.000 

0.6000 

0.0000 

-0.042 + i 1 .89 

-0.057 + i 3.80 

-0 100 


5.6 

-0.9T6 

-o.tfit 

1 , 

0.4028 

-0.080 + i 0.89 

-0.050 + i 1 .86 

-0.037 



3.115 

0.0004 

-0.0115 

-0.036+ i 1.99 

-0.059 + i 3.91 

-0.108 

PARAHETER8: % 

» 3 # 0 1 h 

A- 

* 0 . 1 1 H 

« 0.1. ^ 

: « 0.1, <r» 0.1. 
p 
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•long with th» »olution« in T«bl» 4.t. Tht roll and yaw 
raapon*** obtainad for casas A and B ara dapictad In 
Fig. A. a. Idantical initial conditiona ara uaad for both the 
case* and ara asantionad in th# figure. 

The roll rata in case A build* up saoothly to the design 
value while the roll rat* in case B show* the transient 
r#*onanc* phanoaanon. Thus, th* roll response is as 
expected. However, th* predicted aaplif ication in th* yaw 
response i» hot seen to occur. On the contrary, the 
oscillation* for case B in Fig. a. 4 show distinctly increased 
dasping during the itae th* roll rat* is passing through the 
resonance region. This can be explained by a coaparision of 
the danping ratios calculated froa the complex eigenvalue 
pair* in Table 4.1 which showed th* solution at resonance to 
be more daaped in yaw than the design solution. 

To ensure that th* result indicated by Flg.4.4 was not a 
special case, »# consider another set of paraseter* as 
follow* 

% - 3.0, h - 0.1, - B70*. 

Css* As g * 0 Case B* «.• 5.0 
Again, case A with C equal to tero will have only one 
solution (design solution!. Th# stable solution* in case B 
are listed in Table 4.1 along with their eigenvalues. 
Fig. 4. 5 show# the roll and yaw response* for both case# A 
and B. The roll rate for case A builds up to th# design 
value while for case » It spends •«>«• time in the vicinity 
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of revf-rst resonance. As in Fig. 4. 4, one can observe an 
increase in the damping of the yaw response for case B in 
Fig. 4. 5. The eigenvalues in Table 4.1 for the solutions 
corresponding to the above set of parameters show the 
damping of the resonance solution to be more than that of 
the solution at the design value. Thus, it is likely that a 
response with increased damping would occur for other cases 
of transient resonance too. 

This result is contrary to expectation and has not been 
reported in the literature so far. However, the actual 
response may be complicated by the presence of 
nonlinearities in the pitch-yaw equation (3.6) which have 
been neglected in the present analysis and can be verified 
only from wind-tunnel and flight tests. Induced side forces 
and moments arising in the flight of rolling missiles were 
predicted by Nicolaides C93 and have been described in the 
llteratureC£0-BE3 . These are known to be highly nonlinear 
and could be significant enough to affect the flight of a 
projectile at resonance. Murphy 153 explained how these 
induced side moments could cause the "catastrophic yaw- 
postulated by Nicolaides E93 . 

4.3 Catastrophic Yaw 

Finned projectiles experiencing high-amplitude lunar yawing 
motion are known to suffer from dynamic instabi 1 ity C43 . This 
results in an unlimited build-up of yaw called catastrophic 
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yaw. Many flight failures have been attributedE53 to this 
phenomenon . 

Chadwick Cl 1 3 describes the build-up to catastrophic yaw of a 
missile. The first step is attaining a stable rolling motion 
with The roll rate locked in at resonance. The roll rate 
then shows small oscillations about the equilibrium roll 
rate at lock-in. The lunar motion gives rise to induced 
side forces and moments, which destabilize the yawing motion 
of the projectile. If these induced moments can alter the 
stability in roll of the equilibrium solution at resonance. 


the projectile may 

break out of 

lock-in 

before the yaw 

bu i 1 d-up 

reaches 

a large 

value. 

, This is 

the phenomenon of 

transient 

roll 

lock-in 

and 

break-out 

described by 


Barbara Cl S3 . 

The present study explores the likelihood of occurrence of 
catastrophic yaw even in the absence of induced moments. To 
this end, we first investigate the effect of parameter 
variations on the stability of the equilibrium solution at 
resonance. The base values of the parameters assumed for the 
study are as follows: 

is - 3.0, h » 0.1, fpj « 90 , G « 5.0. 

H « 0 . 1 , Kj,« 0 . 1 , o- ■» 0 . 1 . 

The equilibrium solutions for this set of parameters are 
listed in Table 4.e. Only two of the five equilibrium states 
I corresponding to ll^^l.Oll and ^^=£.861) are stable and the 
eigenvalues for these are also indicated in Table 4 , E ^ ^ 



-TABLE A.1 EFFECT fiF PARAHETER VARIATIONS flN 

EQUILIBRIUM SOLUTIONS 


( 

k ; 

i 

Re{ 

ImC '$5.) 

EIGENVALUES 




H » 0.1, 0 . 1 

H « 0.08, K^- 0 .£ 

-1 . 086 

0 . £469 

0.4079 

UNSTABLE 

UNSTABLE 

-1 . 0E7 

0.7583 

0 . 4041 

UNSTABLE 

UNSTABLE 

1.011 

-0.9439 

0.£045 

-0.087 + i 0.95 

-0.040 + i 1.87 

-0.044 

-0. 178 + i 1.35 

+0. 0£0 + i 1.87 

-0.044 

1 .EA1 

-0 . 0404 

0. 1759 

UNSTABLE 

UNSTABLE 

e.86i 

-0 . 0006 

0.0139 

-0.049 ± i 1 .75 

-0.055 + i 3.67 

-0.088 

i 

-0.075 + i 1.75 

-0. 015 + i 3.67 

-0.177 


PARAMETERS: « 3.0 


0.1, » 90, 6 = 5.0, <r - 0.1. 
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A A 

Parameters H and were selected for the parametric 

A 

analysis. H represents the damping in yaw of the projectile 

A 

and stands for the roll damping effect. The root locus 
plot for the eigenvalues of the solution at resonance 
<^^=1.011) with variation of parameters and H is 

respectively presented in Fig. 4. 6(a) and Fig . 4 . 6 ( b ) . The 
eigenvalues (poles) of the equilibrium solution at resonance 
for the base parameter values assumed above are shown by the 
symbol 'X'. The two pairs of complex poles represent the 
pitch-yaw motion, while the lone real pole represents the 
rolling motion. The arrows show the direction in which the 
poles shift as the parameter value is varied through its 
base value. Only qualitative shift in the pole location is 
shown; a longer arrow indicating a larger shift. 

As is well known, stability requires all the poles to lie in 
the left half plane. Fig. 4. 6 reveals that the complex pair 
of poles closer to the imaginary axis crosses over to the 
right half plane with 

A 

(i) an increase in the value of K^, and 

A 

(ii) a decrease in the value of H. 

This results in the equilibrium solutions at resonance 
becoming unstable. 

To illustrate, we consider the following values of the 
parameters H and Kj, with all the other parameters fixed at 
the samo values assigned to them earlier (Table 4.2): 

H • 0.08, • o.eo. 
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The eigenvalues of the solution at resonance (^g,= 1.0ll) for 
this set of parameters are also listed in Table 4.£. A pair 
of complex eigenvalues can be seen to have positive real 
part indicating instability. On the other handj the design 
solution remains stable while the other three unstable 
solutions given in Table 4.S were seen to remain unstable 
for the new parameter values. Thus, the system has only one 
stable solution, viz. the design solution. Under these 
circumstances, it seems logical to assume that the 
projectile will always end up at the design solution for all 
initial conditions. 

To verify this premise, the response from a numerical 

simulation for this set of parameters is studied and the 

results are shown in Fig. 4. 7. Contrary to our expectations, 

the roll rate can be seen to get locked in at resonance and 

it then shows small oscillations about the resonance roll 
* 

rate <(^*1). The angle of attack magnitude shows a gradual 
build-up. Soon, large amplitude yawing motion is seen which 
lead to instability. This response is typical of 

catastrophic yaw. It may be noted that in contrast to the 
bounded oscillations in yaw seen for cases of stable roll 
resonance (Fig.4.1 and Flg.4.E), the oscillations in Fig. 4.7 
for an unstable resonance solution are unbounded. 

The above analysis reveals two important results. 

(l).It is possible to account for the phenomenon of 
catastrophic yaw without considering the Induced side 





forces end moments in the analysis. This result has not 
been reported in the literature so far. This may serve 
to explain many flight failures where the missile 
unexpectedly destroyed itself after getting locked in at 
resonance despite being designed to avoid roll 
resonance . 

(iil.An unstable equilibrium solution at resonance does 
not rule out the occurrence of lock-in. This directly 
refutes the assertion made by MurphyClSl that unstable 
equilibrium points have no engineering significance. In 
fact, an unstable lock-in solution leads directly to 
catastrophic yaw, whereas for a stable lock-in solution 
the occurrence of catastrophic yaw needs the presence of 
induced moments. In that sense, an unstable lock-in 
solution is potentially more "dangerous". 

4.4 Roll Break-Out 

From the results of the parametric study in Fig. 4. 6, the 
single real root representing the rolling motion does not 
show any change in its location with variation of the 
parameters. As such, the occurrence of break-out from 
lock-in does not seem possible. This was also confirmed by 


studying 

the 

response 

for 

varying 

sets 

of parameters. 

However , 

by 

including 

the 

induced 

roll 

moment terms 

described 

by 

Murphy Cl 93 

in 

addition 

to the 

rolling moment 


due to an offset center-of-mass, it may be possible to 
observe the break-out phenomenon in finned projectiles. 



Chapter V 


CONCLUSIONS 


5.1 Summary and Significance of Results 

(a) Resonant lock-in of finned projectiles is seen to occur 
as a result of aerodynamic and mass asymmetries 
represented by the parameters h, and S. These arise 

from manufacturing and assembly tolerances* fin 
adjustment, and measurement inaccuracies Cl 01 . It is seen 
that lock-in solutions are possible only for certain 
combinations of these parameters. However, these 
asymmetries are essentially random in nature and the 
designer has little control over them. He can, at best, 
ensure stricter control over inaccuracies arising in the 
production process to minimize the asymmetries. This 
will, in general, make a lock-in solution less likely. 

<b) Between normal and reverse lock-in solutions, it is seen 
that larger values of mass asymmetry are needed for a 
reverse solution to exist. This may account for the lack 
of any reported flight experience of reverse lock-in. 

(c) The occurrence of resonant lock— in is largely influenced 
by the initial conditions at launch. These arise as 
result of disturbances during the process of firing a 
projectile and, normally, take the form of a yaw rate. 





For parameter combinations that indicate the existence 
of stable resonant solutions, lock-in for a smaller 
value of 6 was seen to require larger initial 
perturbations. Thus, limiting disturbances at launch in 
conjunction with small asymmetries could be a possible 
solution to avoid lock-in. 

(d) Sixteen sets of initial conditions were investigated for 

each of four different parameter combinations to study 
the influence of the initial conditions on the steady- 
state response of the projectile. The results revealed 
that stable lock-in solutions do not occur frequently. 
It is seen that an unstable solution always exists near 
every stable resonant (lock-in) solution which 

effectively limits the number of trajectories of the 
system leading to the lock-in solution. 

(e) The phenomenon of transient resonance is shown where a 
projectile does not get locked in but stays at resonance 
for some time as its roll rate builds up to the design 
value. Our analysis suggests that projectiles 
experiencing transient resonance show increased damping 
in their yaw response. This result has not been reported 
in the literature so far. However, it is felt that the 
actual behavior of the projectile in flight may be 
complicated by the presence of nonlinear induced side 


forces and moments. 


(f) The present study reveals the significance of unstable 
resonance solutions which had been dismissed as 
unimportant by MurphyC193. It is seen that an unstable 
lock-in solution does not rule out the occurrence of 
roll lock-in. 

(g) Further, a projectile which gets locked in at an 
unstable lock-in solution is seen to exhibit the 
phenomenon of catastrophic yaw. Thus, our study reveals 
that the occurrence of catastrophic yaw can be explained 
without introducing the induced side moments in the 
analysis. 

(h) The phenomenon of roll break-out of finned projectiles 
cannot be observed on the basis of the present study. 

5.£ Recommendations for Further Work 

(i) The analysis carried out in the present study can be 
expanded to incorporate the induced side forces and 
moments described in the 1 iteratureCS0-2£3 . Such an 
extension can be used to investigate the phenomena of 
transient resonance and catastrophic yaw. 

(ii) The possibility of break— out from lock— in of finned 
projectiles can be studied by Including the induced roll 
moment terms described by MurphyC193 in addition to the 
rolling moment due to an offset center-of-aass . 
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APPENDIX t 


Equations (£.1) can be combined to give 

Fy + iF^^ = m t (V + iw) - iu(q + ir) 3 ...(A.t) 

On expressing F^ and F^ in terms of non-dimensional force 
coefficients 

Cy + iCj. « (2m/pu^S) Cf- i(q+ir) 3 - { A . 2 ) 

where it is assumed that V u . 

The linear aerodynamic force for a slightly asymmetric 
missile can be written asC193 

Cy + i * — C ~ expCi(^ + 0^^ ) 3 . . . ( A . 3 } 

(A.S) and tA.3) can be solved to give the complex angular 
velocity as 

~ A, 

q + ir « i ( puS/Em ) C-CjJl - C^expCi(0 + 0^)3 - 

. . . {A.4) 

Equations {£.3) can be put together as 

M+iN = l(q+if)-iplK(q + it') ...(A.5) 

Using non-dimensional moment coefficients 

C^+ (El/pii^SlJ C(q + ir) - iplTiq + ir)3 

...<A.6) 

A linear expansion for the aerodynamic moment can be written 
as Cl 93 

c^+ ic^» 

-iC^^expCi(0 + 0^)3 ...(A. 7) 

where * denotes differentiation with respect to 


dimensionless time, s. 


Using (A. 4) and its derivative in (A.fii along with (A. 7} 
gives the linear pitch-yaw equation asC19!l 
1^" + (H - icrOMp - (M + i<re>'T)'| 

= expEi(a> + ...(A.8) 

where 

Hy^= -(psi^/£l) CC^-ikJ (i-trt®‘c^expci((p^- #^>33 

... (A. 9} 

Equation (3.3) for the roll soment in conjunction with (£.£) 
gives 

LpP + p ..(A. 10) 

Using the expression for L© from (3.4), this can be put in 
the notation of RefI193 as 

0^ ' + C0' - 4/ - iKq (| - 1^)3 » 0 . . (A. It ) 

where ' refers to derivative w.r.t. dimensionless time, s. 

HurphyC193 introduces another dimensionless time,'?', where 
r « E -H/( 1-tr)3'^s ..(A.1£) 

The term within the brackets in (A. IE) can be seen to be 
equal to the resonance roll rate given by (£.17). Thus, the 
effect of the transformation (A. 12) is to scale the roll 
rate in (A.ft) by its value at resonance. Now, the roll rate 
at resonance will be given by . 

With the simplifications adopted by RefC193, (A. 8) can be 


transformed to 
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^ + (H-icr0)^ + ( 1 -<r-iO'0T )$ « 1 -cr)exp Ci ((J+<|^3 

..(A. 13) 

where denotes derivative w.r.t. the independent variable 

T. i in (A. 13) stands for the magnitude of the initial trim 

T-O 

and or is the ratio of axial to transverse moment of inertia. 


At the steady-state, the nutation and precession modes are 
aesumed£193 to have damped out and the response is taken to 
consist only of the trim. 

I = fe,.exp( i<t)) 

where substituting ( A . 1 -4 ) in (A. 13) gives 


(A. 14) 


^ ■ TTlFTTTr 


)] 


(A. 15) 


At resonance (I^I*!), the trim angle magnitude can be 
written from (A. 15) as 

. . (A. 16) 


Using (A. 16), a scaled complex angle of attack can be 
defined as 

^ » I ..(A. 17) 

(A. 11) and (A. 13) can be put in the following form using 

^ ci - lei's -¥>3*0 . .(A. 18) 

^ + CH + iC£-0l^35 + C 1 -er) Cl -♦)^+i<|>h + 1^*35 

“ ( 1 -or) f h jexpl i<^^^) ..(A. 19) 
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